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SUMMARY 


The work reported was undertaken in connection with 
a vibration study of helicopter rotors. The analysis of 
vibration and certain other problems require a knowledge 
not only of the average induced, velocity but also of its 
distribution around the rotor disk, A concept of a 
simplified vortex system of a rotor wake is used in 
obtaining a. formula for the normal component of induced 
velocity along the fore-and-aft diameter of the rotor 
disk. This formula is intended to represent the main 
effect of the skewed, wake in producing an uneven velocity 
pattern at the rotor. This induced velocity is expressed 
in terms of elliptic integrals as a function of the skew 
angle and of the circulation per unit axial length of the 
wake. A simple approximation to this function consisting 
of the first two terms of a Taylor expansion and giving 
the value and the slope of the induced-velocity function 
at the center of the disk is also presented. 

An approximate method of representing the induced 
velocity in terms of flight velocity for a given thrust 
has also been indicated by combining the present theory 
with Crlauert’s thrust equation. 

A comparison of the theory with corresponding values 
of induced velocities computed by Seibel's formula from 
p i t c hi. ng -moment data pub 1 i s he d by Whe a t le y and Bi o 1 e 1 1 i 
is taken as evidence that the most significant factors 
have been taken into account. 
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INTRODUCTION 


The work reported herein was undertaken in connec- 
tion with a vibration study of helicopter rotors. Although 
for certain purposes the use of an average induced velocity 
assumed constant over the rotor disk gives reasonable 
accuracy, the analysis of vibration and certain other 
problems require a knowledge of the distribution of Induced 
velocity around the rotor disk. Very little such informa- 
tion, however, is available. Several writers have recog- 
nized that the induced velocity increases toward the rear 
of the rotor disk but have left the magnitude of this 
increase as a matter of speculation. Glauert (reference 1) 
and Wheatley (reference 2) have used an induced-velocity 
formula containing a first harmonic term with an unknown 
coefficient K having a value somewhere between 0 and h 
Nikolsky (reference 3) has suggested a formula that is 
equivalent to putting K - 2 \i, where p is the ratio 
of the flight velocity parallel to the disk to the tip 
speed. Seibel (reference ).j.) has recently suggested a 
method of deducing values of the first harmonic induced 
velocity from experimental measurements of the pitching 
moment of a gyroplane by Wheatley and Bioletti (reference 5 ) • 
Seibel gives an example of the use of such data for 
assumed values of pertinent parameters and shows that the 
resulting variations of induced drag may be expected to 
produce a peak of vibration at around 25 miles per hour 
for a typical case. 

In the present paper a concept of a simplified vortex 
system of a rotor is used for obtaining a formula for the 
normal component of induced velocity along the fore-and- 
aft diameter of the rotor disk. 


BASIS OF ANALYSIS 


The aerodynamic theory of helicopters can be thought 
of as a generalization of propeller theory, thus certain 
analytical methods found in propeller literature can be 
extended to apply to helicopters. The method used herein 
is based on the assumption of an idealized picture of a 
wake vortex pattern consisting of an elliptic cylinder 
which Is, in general, skewed with respect to the rotor 
axis at an angle that depends upon the flight velocity 
and upon the Induced velocities. Enough simplifying 
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assumptions are made about the wake vorticity distribu- 
tion to obtain a tractable form for the integration of 
the Biot and Savart law and thus to represent the main 
effect of the skewed wake in producing an uneven induced- 
velocity pattern at the rotor. 

The wake pattern is assumed, for the present purpose, 
to form a continuous distribution of vortex lines on the 
surface of an elliptic cylinder such as would be formed 
by the vorticity shed from an infinite number of blades 
with constant circulation and very light loading. The 
analysis is further simplified by making use of the 
property that the induced-velocity field of the helical 
vortex system can be considered formed from two simpler 
vortex systems of which one is composed of circular vortex 
rings and the other of axial vortex lines. (See fig. 1.) 
The circular rings determine the fore-and-aft distribution 
of velocities, which are the important ones at small 
advance ratios; whereas the axial lines determine the 
rotational velocities of the wake, which are not con- 
sidered important for the present problem. 

With the assumed picture of the vortex wake, the 
normal component of induced velocity along the fore-and- 
aft diameter of the rotor disk is obtained in closed 
form, as a function of the wake geometry and of the 
strength of the vortex sheet. 

Although the skew angle and the vortex strength are 
the fundamental variables of the present problem, it is 
desirable in applying the theory to have the induced 
velocity expressed in terms of flight velocity for con- 
stant thrust. One method of representing the results in 
this form has therefore been indicated by combining them 
with Glauert’s thrust equation. 


SYMBOLS 


% 


w ake s ke w an g 1 e 


5 , 


X, 


V 


z 


p, t coordinates relative to 
y, z coordinates relative to 
n o rm a 1 in due ed-velocity 


s ke w e d w a ke (see fig. 2 ) 
rotor disk (see fig. 2) 
component at rotor disk 
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v, ( ' , v * , v £ », v t » Induced- velocity components in 

^ ultimate wake 

u value of induced velocity at center of rotor disk 

u-i rate of change of induced velocity along fore-and- 
aft diameter of rotor at r =0 

u hovering induced velocity 

V flight velocity 

V f resultant velocity at rotor disk in Glauert ? s 

formula 

a angle of flight path to rotor disk 

r distance from huh to position on blade; expressed 

nondimensi onally in terms of rotor radius 

v|/ azimuth, angle of blade position 

T thrust 

p air density 

A area of rotor disk 

C T thrust coefficient 

dZ line element 

ANALYSIS 


Consider the vortex wake pattern of figure 1(b). 
Appendix A shows that, if the vortex lines are circles, 
the induced velocities inside the ultimate wake (see 
fig. 2) are related by the equation 




( 1 ) 



tt - % must- be used in place of • 
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The evaluation of the induced velocity at the disk 
is given in appendix B. The complete expression for the 
distribution of the component normal to the disk along 
the fore - and -aft d i ame ter is 



~ K(r) 

tt sin ^ 


r(r + 1) 


/ p p 

l 4 .Tr tan^ \/tan L p^+ r 




where 


( 2 ) 


K ( r ) and 





are complete elliptic integrals of 


the first and third kinds , respectively, and , b^ , 
and i are defined following equation ( B12 ) in appendix B. 

For many purposes a sufficient approximation to 
this expression is obtained by using the value and the 
slope of the function at r = 0 representing the first 
two terms of a Taylor expansion 


v z = u + u^r cos f + 



where, for % < 


and, for 


% > 2 ’ 


U- 


u tan 


% 

2 


u 


1 


u 


cot 


% 

2 


That the higher-order terms are small is indicated, 
from, the plot of the exact expression for vertical induced 
velocity against radius for the extreme case in 

which % - g-. (See fig. 3 .) 
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In order to make use of the expressions developed 
herein, the skew angle 9 C must he related to other flight 
parameters. A simple method is to combine equation (3) 
with the expressions given by Wald (reference 6), which 
are based on Glauert’s formula for thrust. The Induced 
velocity u in Glauert’s formula 


T = 2 PAuV T (]+) 


Is Identified with the value at the center of the disk 
in the present theory. 

From the geometry of the wake, it can be seen that 


T 

c 

~ V 


cos 4- a) 


( 5 ) 


Then by using the relations. 


for % < § > 


v £* 

Y 


tan 


% 

2 


( 6a) 


for 



and 


it follows that 


v 



cot 




> 


J 


(6b) 


( 7 ) 


% 

tan 2 


V 

u 


cos(^ 4- a) 


( 8 ) 
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I f Wald. f s e qua 1 1 on 


u^ + 2 My sin a + 


" u 


o 


(9) 


is written in the form 


/hi 

V a c 


k 


y 


i 


1 + 2 — sin a + (~ 

u V u 


2 


_V_ 

u„ 


_irV 
u u 


( 10 ) 




equations ( 3 ), (3) ? and ( 10 ) can be 
of u/u Q and u^/uq against v/u 
These plots are shown in figure ii. 


used to obtain 
for constant 


plots 

a. 


These curves indicate 


that the first harmonic velocity term is a maximum at a 
fairly 1 ow f 1 i ght speed, as w as expected from Se i be 1 f s 
interpretation of pitching-moment data. The curves also 
show that the height of the peak decreases as the angle a 


increases toward 90 °, 
of a helicopter. 


corresponds to vertical climb 


values of 

. 4- 


Curves are also plotted in figure I 4 . for negative 


a but, si 


ice the vortex ring 


and the turbulent 

states of operation are included in this range, the vortex 
theory based on the assumption of a long cylindrical wake 
cannot be expected to give reliable results. Part of the 
curves of figure [ 1 . have consequently been plotted as 
dashed lines to indicate their provisional" status. 


A direct experimental check of the aoplicabili 
the present theory has not been possible because of 
unavailability of data on induced velocities. A te 
comparison with Seibel ? s method has been made, howe 
by computing the induced velocities corresponding t 
series of test conditions for which pi tching -moment 
are available and then using Seibel’s formula to co 
values of induced velocity for comparison with the 
of the present theory. 


ty of 
the 

ntati ve 
ver , 
o a 
d.ata 
mpute 
results 
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A plot of u i/ u 0 against V / u Q , determined by 

Seibel’s formula and the data of reference 5? is shown 
in figure 5 for two values of the blade pitch setting. 

For comparison, the formula for u^/u 0 given in the 
present paper has been applied ( by use of the 

formula — = — J to the same test conditions 
u cos a | Cm/ 

o v V 

specified by the values of a, C T > and p in refer- 
ence 5* The resulting curves are also plotted in fig- 
ure 5? which shows that both methods give peaks of the 
same general appearance. Systematic errors influencing 
the magnitude of the theoretical peak may be expected 
from several sources, which include assumptions of 

(1) Vortices shed only at blade tips 

(2) No contraction of wake 
(5) Infinite number of blades 

(]q) Slope of induced- velocity function as repre- 
sentative of front and back inequality of 
veloci ty 

The net effect of these assumptions appears to be to 
underestimate the magnitude of u-,. The use of an 

effective radius smaller than the actual rotor radius 
would give larger induced velocities for the same total 
thrust and would tend to show closer agreement between 
the curves. 

In view of the various approximations and possible 
sources of error involved in both methods, the qualitative 
agreement between these two methods is evidence that the 
most significant factors have been taken into account. 
Further testing and refinement of analysis should lead 
to a more detailed understanding of the phenomena 
involved. 


CONCLUSIONS 


From theoretical considerations. It was concluded 
that the Important variable determining the increase of 
induced velocity toward the rear of a helicopter rotor 
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in forward flight is the angle ^ by which the wake 

axis is skewed from the rotor disk axis. If the induced 
velocity distribution is represented by a series 
expression of the form 


v z = u + u-^r cos \| / + ... 


integration of the Blot and. Savart law for a. simplified 
form of vortex wake gives the value 


u^ - u tan -p- 


where 



v z 

normal component of 

i n due e d ve 1 o c i t y 

u 

value of v at center of disk 

r 

nondimen si on al r-adia 

1 coordinate 


azimuth angle 



The exact expression for the normal component of 
the induced velocity along the fore-and-aft diameter for 
the assumed vortex wake was also found in closed form 
in terms of complete elliptic integrals of the first and 
third kinds. 


It i 
use of an 
using an 
radius . 


S o believed that systematic errors due to the 
idealized wake may be partly compensated by- 
effective radius smaller than the" actual rotor 


Langley Memorial Aeronautical Laboratory 

National Advisory Comm.it tee for Aeronautics 
Lan g 1 e y Pi e 1 d « V a . 
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APPENDIX A 


VELOCITIES IN TEE ULTIMATE WAKE 


The assumption of circular vortex lines on the 
surface of an elliptic cylinder forming the ultimate wake 
is shown to imply a velocity which is constant inside 
the wake and. such that 


vj 

_JL 

V 


tan 


X 

2 


The method Is first to try to satisfy all the require- 
ments of the problem by assuming that inside the cylinder 
the velocity components v^h and v^* are constant but 

of undetermined magnitude. Their constancy means that 
the wake cylinder moves as though it were a rigid cylinder 
and the flow outside is consequently the same as the 
well-known velocity field outside a solid cylinder. 

The vortex strength can then be determined in terms of 
the assumed, velocities by evaluating the line integral 
of velocity around a suitable path enclosing a portion 
of the surface of the cylinder and equating this Integral 
to the total circulation of vortex lines that thread the 
path of Integration. This equation yields an expression 
for the slope of vortex lines on the cylinder, which can 
be Integrated, to give an equation of the vortex lines. 
Comparison with the equation for circles shows that 


Vpe* 

V 


tan 


D 

2 


Consider the line integral of velocity around, path A 
of figure 6. If the circulation per unit axial distance 
along the wake is denoted by P/s, the line integral 
around path A shows that 




t 
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N o w go n s i & er the v e 1 o c i ty i n t s g r a 1 
the gr -plane of figure 6, Referred to 
air 5 the velocity potential on the outs 
motions in the Bp -plane Is 


for path B in 
axes fixed In still 
de surface for the 


0 = -v ? a cos 0 

where v r ' refers to the constant velocity inside the 
e, 

cylinder* This potential Is the same as in the flow at 
the surface of an elliptic cylinder moving through an 
incompressible fluid, that is at rest at Infinity. Then 

/v . d Z — gp- d0 - v*h d£ T 

U OU b 

= v^ ? (a + b) sin 0 dQ 

in 

This Integral is equal to the total circulation threaded, 
at B. If the slope of the vortex lines Is represented 
by d£/de, the number of lines that thread the path B 
i s 


i ciL 
s dG 


19 


Consequer fcly 


v , ? ! { a + b) sin 9 


and * af t e r I nte g r a t 1 on . t h 1 s e a ua 1 1 o n b e c o me s 


V V ' ti 


( s + b ) c o s 9 


The equation of a circular vortex, line , however* is 


£ ~ -c cos 0 
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where 




Henoe , 


v- * s - v^> ? c 

3? V <■ * fiL 4 * h 


sin % _ 
1 + cos ^ 


= tan - 


% 


2 


It should, be noted that this relation implies that, 
insi de the uiti mat e -wake cyli ride r , 
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APPENDIX B 


INDUCED- VELOCITY .FIELD AT ROTOR DISK 


As stated previously, the entire vortex wake is 
considered to be an elliptic cylinder composed of a 
continuous distribution of circular vortex filaments 
parallel to the rotor disk. (See fig. 7 • ) By using 
the Biot and Savart law, the velocity vector induced by 
the vortex wake at any point P having coordinates 
(r,^ ) in the plane of the rotor disk Is given by 


v = 



i 3 k 

r cos\J/+cos 8 -z tanX -sin G ~r sin\j/ s 
sin 6 cos 0 0 


dz d9 


where 

P/s axial-vortex strength 

r nondxmensional distance in terms., of rotor radius 

m = tan "X 

h k unit vectors along x, y, and z directions, 
respectively 

p distance from point (r, ) to element d 0 dz of 

wake surface 

The value of p is given by 

p2 5= 2 + r^+ z^ (1+ m2 ) - 2rzm cos -2m cos 0 +2r cos(7 - 0 ) 


p' 

(Bl) 


The vertical component of the induced velocity v 7 at 
the point (r,\[r) can be written from equation (Bl) as J 


i4 
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1 - zm cos 9 + r cos (\|/ - 0) 

— — - - - dz do 

p? 


(B3) 


The value of v z along the fore-and-aft diameter of tfie 

rotor disk can be determined by setting \[/ — tt in 
equati ons ( B 2 ) and ( BJ ) ; thus 


V - V 

z , z 




2tt 


00 


1 +111 
Utr 


[l ~ ( zm+ r) cos 0] dz d9 


0=0 c z ~0 


l + r 2 + z 2 (l + m 2 ) + 2rzm - 2( mz + r) cos 0 


.13/2 


(bU) 


where F/s has been replaced by the vertical induced 
velocity in the ultimate wake v z h (See appendix A.) 

The integration with respect to z can be readily 
performed by use of formula 200 of Peirce’s n A Short 
Table of Integrals.” The result is 



1 + 


r cos 9 - r‘ 




1 + r 2 -*■ 2r cos 9 + sin^0 m \J]_ + - 2r cos 0 


r\/l~+ m^ yl - 2r cos' 9 + 
m (l + - 2r cos 9 + sin^o) 


d0 


(B5) 


in which the limits have been changed to 0 to rr since 
the integrand is an even function of 9. 

The value of r is assumed to be positive in the 
integrations. For negative values of r, the substitu- 
tion 0 = it - 0 converts the integral Into one similar 
to the integral for positive values of r. The Integral 
of the second term of equation (B5) is shown later to be 
identically zero for all values of |r| < 1. It can 
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then be seen that u, which is the average value 
of v z across the fore-and-aft diameter, is equal to 

the value of v z at the center of the rotor disk since 
the curve of v z against r is symmetric about r - 0. 
If r Is put equal to zero in equation (B5), 


( b6) 


The rate of change of v z along the fore-and-aft 

diameter can be immediately found from equation (B5) by 
differentiating within the integral sign (this 
differentiation is permissible since the denominators 
are always greater than zero for |r| < 1), setting 
r = 0, and then evaluating the resultant integrals. The 
result is 
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The ccmplete evaluation of equation ( B5 ) is performed 
most readily by integrating each term separately. The 
integration is restricted to the case of |r| < 1, that 
is, within the rotor disk. 

The integral of the first term of equation ( B5 ) is 
evidently rr. 


It is now shown that the integral of the second 
term of equation (B5) is identically zero. By using the 

n 

conventional transformation tan ~ = z, the integral 
becomes 



( cos 8 - r)d6 - 2r I [1 - z c - r ll + z^Jj dz 

1 + r 2 - 2r cos8+n 2 sin 2 0 t / _ q + r 2) + Z 2A 2 _ 2r A . ^ + ^ 


r 


1- z 1 2 - r(l+ z 2 ) 


2 2 
z 


(B9) 


The integrand is finite everywhere, approaches zero as z 
becomes infinite, and is zero and changes sign only at 


the point 
parts ; 


1 - r 
1 + r 


Divide the integral into two 


oo 



dz 




dz (BIO) 


1 - r 1 

If the transformation z = - is applied to the 

1 + r w 

second integral of the right-hand side of equation (BIO), 
this integral is seen to be identically equal to the 
negative of the first integral. The total integral of 
equation (BlO), or equation ( B9 ) , thus is zero. This 
integral Is not zero, however, if |r| >1. 
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The integral of the third term of equation (B5) is 
equal to 


ryl + iri 
m 


2 




e-o 


de 


1 + r - 2r cos 8 


2 ry 1 + m‘ 
m 


2 


K(r) 


where K(r) is the complete elliptic integral of the 
first kind. 

The integral of the last term of equation (B5) is 
evaluated by first using the substitution tan i. = z, 

then separating into partial fractions, and thus 
obtaining 


L | 2 

yl + m 


.TT 


1 - 2r cos 8 + r 


2 


9=0 1 ” 2r cos 8 +r ^ + tp? sin^8 


.cL8 



2rrym^ + r^ 


where 


2 2 
z + m 2 


k 2 = /. L - Z .-., £ 
vr + r 


2 


dz 


( Bll ) 


m- 


+ 1 + 
r + 1 


m 2 


1 - ym^ + r^ 


r + 1 
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It is noted that > 1 and m, < 1 for a ll positive 

values of r less than 1. By using the substitution 


Q 

z + 1 


the integral of equation (Bll) is trans- 


formed into the following form after again being 
separated into partial fractions: 



b i = 


Integral ( B12 ) is immediately transformed into complete 
elliptic Integrals of the third kind by multiplying the 

numerators and denominators by \/l - j^y 2 and separating 


into elliptic Integrals of the third kind. 


By collecting the various integrals, the complete 
evaluation of v„ for positive values of r is 

h=ir 


-V = K(r) 

v z 2 tt sin ^ 


( b i + f) n d' b i’ j) + ( b 2 - j 2 ) n(|.-t 2 , j) ( B13) 


For negative values of r, the sign of the second 
and third, terms of equation ( B 13 ) should be changed and 
r should be considered a positive number throughout the 
resultant formula. 


r( r + 1) 


i+Tf tan pd\/tan^+ r‘ 
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Figure 1.- Representation of skewed helical vortex wake by circular and linear vortex 

wakes . 
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Figure 5.- Comparison of theoretical indue 

from pitching-moment data 
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